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A toy model of Einstein gravity with a Gauss-Bonnet classically "entropic" term mimicking a 

fT^ , quantum correction is considered. The static black hole solution due to Tomozawa is found and 

generalized with the inclusion of non trivial horizon topology, and its entropy evaluated deriving the 

("^ ' first law by equations of motion. As a result the Bekenstein-Hawking area law turns to be corrected 

04 , by a logarithmic area term. A Misner-Sharp expression for the mass of black hole is found. Within 

» I . a Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmological setting, the model is used in order 

^ l' to derive modified Friedmann equations. Such new equations are shown to reproduce the first law 

.^^ \ with the same formal entropy and quasi local energy of the static case, but here within a FLRW 

space-time interpreted as a dynamical cosmological black hole. A detailed analysis of cosmological 

>0 ' solutions is presented, and it is shown that the presence of the correction term provides regular 

solutions at the origine and interesting phases of acceleration and decelerations, as well as, with 

' fj', negligible matter, exact de Sitter solutions. 
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I. INTRODUCTION 

QQ \ In a recent paper [l|, Tomozawa has put forward an argument to deal with finite one- loop quantum 

r — ■ corrections to Einstein gravity, by reconsidering an old result of himself. The techniques within his proposal 

OO \ refer to are well known and can be found in the celebrated texbook [2] . Here, we present a simplified approach 

motivated by alternative arguments, and in some sense, also inspired by the entropic approach proposed in 
T- H ' Refs. [3, |j]; where Lagrangian surface terms are taken into accout. 

f^ , The paper is organized as follows. In Section II, we will give a short review of the traditional semiclassical 

CO ' approach to quantum corrections of Einstein's gravity. In Section III, the "entropic" approach is introduced 

and the spherically symmetric static case is investigated. In Section IV, the model is extended to the 

spherically symmetric dynamical FLRW space-time and in Section V some explicit solutions are presented. 

Conclusions are given in Section VI. 



II. BRIEF REVIEW OF QUANTUM CORRECTIONS TO GRAVITY 

Before discussing the model which mimics the Tomosawa proposal, we review the more traditional treat- 
ment as proposed by Starobinsky in the seminal paper [6| . The argument is based on the so called semiclas- 
sical gravity approach, where the backreaction of quantum fields are taken into account, in order to correct 
the classical Einstein equation, namely 

Gij = Rij - -gij R = x{Ttj) , X = SttGat , (1) 
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gij, R and Rij being respectively the metric, the scalar curvature and Ricci tensor, Gjv the Newton constant, 
and (Tij) the vacuum expectation value of quantum energy stress tensor, which in general is not explicity 
known. As usual, units measure are chosen in such a way that the speed of light is equal to one. 
By taking the trace of equation above one has 



R = -x(Tl 



(2) 



in which the stress tensor trace appears. When one is dealing with a conformally coupled quantum field, a 
quantum conformal anomaly is present (see the review paper [7|). In four dimension it reads 

IJi)^- {aW + PG + 5AR). (3) 

Here W is the "square" of Weyl tensor Cijrs and G the Gauss-Bonnet topological invariant. They read 



W = c'^^'^a. 



G = R'^f'^R, 



■ijkl 



iR'^R,, + R^ 



(4) 



The coefficients a, /3, and S depend on the number of conformal fields present in the theory. In some 
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-^, N being a very large parameter. 



conformal field theories [8J , one has S = and a = 

One might try to solve the anomaly driven trace equation ^ within the static spherically symmetric 
(SSS) Ansatz where metric has the form 



ds' 



-a{r)df 



a{r) 



r^ da^ . 



(5) 



The account of conformal anomaly contribution to spherically-symmetric space-time solutions (SdS BH or 
wormholes) has been done in Ref. [Ifll (see also Ref. [ll| for the case of extended theories of gravity). In 
such a case one has 



1 d^ 



1) 



G = 



2 d^ 



a- 



1], 



W = 



j,2 ^j,2 L J ) 2 

As a result, the quantum corrected Einstein trace equation reduces to 



dr'^ 



a-1 



dr"^ 



[r2(a-l) + 2x/3(a-l)2] + 



Xa 



dr^ 



a- 1 



0, 



(6) 



(7) 



This is a second order non linear equation in the unknown variable a{r). Its exact solution appears to be 
difficult to find. However, one can directly verify that there exists the de Sitter solution 



a{r) = 1 - 



2x/3' 



(8) 



which correponds to the famous solution found by Starobinsky in FLRW coordinates. 

On the other hand, putting S, = XQ^ = X/^i one can find perturbative solutions assuming ^ to be a small 
quantity, that is 



a{r) = ao{r) + iai{r) + i'^ a2{r) + ... 
Starting from ^ = one finds the unperturbed solution 



Ij^^'-'H 







a{r) ^ ao(r) 



1 _ £l _ £H. 
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(9) 



(10) 



ci,C2 being constants of integration. If we choose ci — 2M and C2 — 0, then such a solution may be 
interpreted as the Schwarzschild solution generated by a body of mass equal to M. 
Now, by taking into account of first order perturbation we get 



a{r) r^ ao(r)+^ai(r) = 1 - 



Cl + <^C3 C2 + CC4 4^C? 8^ClC2 



22ec| 

5^6 



(11) 



C3,C4 being further constants of integration. In principle, one may investigate this quantum corrected 
approximate static solution, but this will be not done here. 

In Ref. [9| (see also [l2|), an attempt has been made to solve the anomaly driven trace equation, but in 
the special case of negligible a. In fact, for a = and £, — k/3 one obtains 



«W = i-^ 



i±Wi + ^ (^ + 1 + "'" 



(12) 



where Ci,C2 are arbirary constants. From the mathematical point of view this is a quite interesting result, 
because it might provide quantum corrections to the Schwarzschild solution, but unfortunately it is not 
physically relevant since a and /3 are quantities of the same order. For this reason such an approach does 
not seem completely adequate to describe physical situations. 

III. ENTROPIC CORRECTED STATIC SPHERICALLY SYMMETRIC METRIC 

Let us come back to Tomosawa proposal, reformulating it within a classical Lagrangian approach and 
working first in n dimensions, and then making a "entropic" dimensional reduction to n = 4, also in the 
spirit of other different approaches (see, for example, [Sl-Q)- 

It is well known that in a four dimensional manifold the Gauss-Bonnet term does not contribute to the 
classical equations of motion. At quantum level, the situation changes due to the regularization procedure, 
and this is the key observation made in Ref. [l| . In order to mimick quantum corrections due to Gauss-Bonnet 
invariant or alternatively in order to activate such a term in the variational principle, we will make use of the 
fact that the functional variation of the classical action in arbitrary n dimensions does not commute with 
the limit n -> 4. 

Thus, let us consider the following classical action 

^(i^-2A+^), (13) 

where A is the cosmological constant, G the Gauss-Bonnet invariant in n dimensions and ^ an arbitrary 
parameter, which will be assumed to be proportional to x/3. We shall derive the field equations and at the 
end of calculation we shall perform the n — >■ 4 limit. 

We shall look for static, spherically symmetric solutions (SSS) with arbitrary horizon topologies. The 
generic metric reads 

ds^ = ^a{r)b\r)dt^ + -— + r^ dJ^l , (14) 

a{r) 

where a, b are arbitrary functions to be determined and dS^. is the metric of a maximally symmetric n — 2 
dimensional manifold, which will reduce to the metric da^ of a two dimensional maximally symmetric 
manifold in the limit n — > 4, that is 

dT,l ^ dal = 4-^ + p^d0^ , for n ^ 4 . 

1 — kp-' 

Here, fc = 1, 0, —1 respectively for spherical, toroidal and hyperbolic topology horizons. 

Within this static Ansatz, the Lagrangian becomes a function of a[r], b[r] and their derivatives, that is 

I^n f drr''-^L{a,a',a",b,b',b"), (15) 



where fi is a constant factor due to the integration on all variable except r. 

Now, by means of the Weyl's method discussed for example in Ref. [lj|, one obtains the field equations, 
which in the limit n — > 4 are regular and read 



1 d 

r^b dr 



{k-a)r-—- + ^ '- 

6 r 



= 0. (17) 



These equations can be obtained by an effective classical Lagrangian, which mimics quantum corrections, 
obtained as the hmit n — > 4 of the one in (|15l) (apart total derivatives). 

Equation ([16]) has the trivial solution b ~ Const, and we can choose 6=1 without loosing generality, 
while equation pT)) gives 




a{r)=k + -\l±^ll + A^[- + -]] , (18) 

where C is a constant of integration. Equation (|18p represents the topological generalization of Tomozawa 
black hole solution discussed in 1] in the presence of a non vanishing cosmological constant. It is the anologue 
of black hole solution of Lovelock gravity in five dimensions [iJ, Il3- Furthermore, it is formally identical to 
the black hole solution found in [la] , as a particular limit of Horawa-Lifsits gravity considered in [iTJ • 

What about the meaning of integration constant C? A possible approach consists in discussing the limit 
5 — ^ of solution (|18l) as in \V\. The limit is finite only for the solution in ([18)) with the minus sign in front 
of the square-root. In such a case the result is 

a{r) ^k---^-. (19) 

r A 

A a consequence one may conclude that C = 2M, M being the mass of black hole. 

Another approach consists in the investigation of the Clausius relation dM — T^dSn, which relates mass, 
temperature and entropy associated with the given black hole solution, as a direct consequence of field 
equations. As we shall see, this approach has the merit to give informations about the entropy Sh of the 
black hole solution. 

First let us study under which conditions the solution (|18|) represents a black hole. As it is well known, 
one should have real positive solution of a{ri{) — 0, namely 

^P - CrH + krl ^^r% = Q. (20) 

The general solution of this quartic algebraic equation is quite complicated and we will not write down 
its explicit form. For our purposes it is sufficient to know that, depending on parameters, there exist real 
positive roots. For example this is certainly true in an asymptotically flat manifold, that is for A = and 
k = 1. In such a case one gets 

(21) 

which is positive for ^ < C^/4 and so it defines an event horizon. 

Associated with any event horizon there exists a Hawking temperature given by 

Th-^, >^H=°f =^ a'H^^nTH, (22) 

Kh being the surface gra vity related to the horizon at rn ■ This result is robust and can be derived by several 
alternative methods [l8l - [2q |. 

Following [21| we evaluate the equation of motion ([T7)) at the horizon, where a// = and a'^ = AttTh- We 
obtain 




AttTh irH + ^]=k^Arj,~^^^, (23) 



2/cA , . , ^k^ dC 
rH ) rjj drn 

where the latter expression is a direct consequence of (j20p . obtained by assuming A to be a given parameter 
and the integration constant C = C{rH) to depend only on the horizon radius. 



Now, introducing the horizon area Ah = Vkrf^, Vk being the measure of the unit surface (for example 
Vi = 47r) we can write the latter equation in the "thermodynamical" form 

It is quite natural to interpret this identity as the Clausius relation TudSu — dEn for the black hole solution 
with entropy Sh and quasi-local energy Eh evaluated on the horizon given respectively by 

SH^^ + 27rk^\n-^, (25) 

and by 

E« = :^ = ^('^^^~r'" + —)- (26) 

2Gn 2Gn \ 3 rn J 

This result gives the Energy of the black hole as classical Misner-Sharp mass plus a correction which depends 
on the parameter ^. Furthermore, with regard to the corrected black hole entropy, the area law {Sh = 
ttAh /VkGisf) does not hold, because a quantum logarithmic correction is present. This is a well known 
general result, and it has been derived many times. It has been proposed in [QJ. and it appears in the 
quantum field theory treatment of black hole entropy with heat kernel techniques [22, l2j| , or loop gravity 
[2J| or other approaches (see the recent paper [23] and the references therein). 

It should be noted that for toroidal black hole, present when A < 0, the correction is absent while for 
spherical and hyperbolic black holes the corrections have opposite sign. 

The expressions (|26|) and p5|) are depending on quantities like vhjAh, and xh which are scalars for a 
generic spherically symmetric (dynamical) space-times. It also follows that Clausius relation should have 
this covariance form. We will verify this fact in a dinamical spherical symmetric space-time as the FLRW 
one. 



IV. ENTROPIC CORRECTED FLRW SPACE-TIME 



Anomaly driven FLRW models have been already considered in the past [6|, |28| - |31| . Here we shall consider 
the toy model described by the action (|28|) . but in n-dimensional, spatially flat FLRW space-time defined 
by means of the metric 

ds2 ^ ^e^mdt' + a{tf{dxl + ...dxl^,) , (27) 

where ri{t),a{t) are arbritary functions, which will play the role of Lagrangian coordinates. The quantity 
'ri{t) will be set equal to zero at the end of calculations. With this choice the parameter t will become the 
standard cosmological time. 

In the total action we must also include classical matter described by a perfect fluid with density p{t) and 
a pressure p = wp{t), < w < i being a constant. Thus the total classical action will read 

/.^/d".y3^(i?-2A-^)+/., (28) 

Im being the action of matter which assumes the form 
1 



2 , d"a;V^.g'^ [{p + p)u,Uj + {p - p)g,,] . (29) 

Making the variations with respect to rj{t) and a{t) in arbitrary dimension n, then taking the limit ri — > 4 
and finally putting 7]{t) = we obtain the two equations 

H^^^ + ^+^H\ (30) 



3 rr2 XP , A , rTT2 ( c^tt . ^ 



H+-H' = -^ + -+^H'[2H+-H'j . (31) 

As usual H — ixja is the Hubble parameter. An equation similar to pop has been obtained in a covariant 
renormalizable model for gravity [27| . 

As a crucial consistent check, it is easy to show taht the latter equations give rise to matter continuity 
law in agreement with difFeomorphism invariance that is 

p = -3i/(l + w)p, (32) 

and also to the generalised Raychoudury equation 

j2 _ X(l + 3u;) p A 2 /^ I- 2 



& = H + H-' = -^^- LL + -+ e^H^ {2H + H^] . (33) 

6 3V/ 

In order to analyze the physical consequences, one may indifferently use two among the three equations 
([5(I)) - (|32p . An equation formally identical to (|30l) has already been obtained within AdS/CFT holographic 
correspondence in 1 35 ..36.] . in brane cosmology ^37.] or assuming Clausius relation and a logarithmic correction 
to the area law in [38|. 

The FLRW admits a dynamical trapping horizon (Hubble horizon) in the sense of Hayward [S^] located 
at Rh — jj, with associated surface gravity and "dynamical temperature" [33J 

«.-,.,^,, r„.|f.-i-(ff + ^J^ (34) 

The identification of temperature with the suface gravity as it happens in the static case, is supported by a 
tunneling computation via the Hamilton- Jacobi method (see, for example, [3j, |3J] and references therein). 
Now we introduce the generalised Misner-Sharp energy evaluated on the dynamical Hubble horizon, by 
means of equation ((26| obtained for the static case, then we get 

1 fr. A , e 




and similarly for the Entropy, 

■'"-M;-"'^'''^- A„^,,B:i.^'^_^ (36) 

Such quantities satisfy the Clausius relation as a consequence of field equations ([30 | - ([33)) . In fact one has 
dEH - ThcLSh^ - ^T^2)dVH . (37) 

where Vh = AirR^/S is the volume of the Hubble sphere, while T(^2) = P ~ P = {w — l)p is the reduced 
stress tensor trace, a scalar quantity in a dynamical symmetric space time [32l |. As for the static case, the 
cosmological constant is considered a given quantity. 

V. EXPLICIT SOLUTIONS 

Here we assume p > 0, A > and analyze the behaviour of the solutions for all possible values of £,. For 
simplicity we also use units for which x = 1. 

First of all, in the absence of matter, namely when p{t) = 0, there exist de Sitter solutions for which 
Hds = 0. In fact, equations (pO| and (|3T|) have the trivial solutions 





1±^/1-^), A>0,C>0, ^CA<1, (38) 



H^ 



''~m[ ' 



4|e|A \ 
3 j' 



A>0,C<0. 



(39) 



In the limit |^A| <C 1 one has (for expanding universe, H > 0) 



i?l~ 


[a 

V 3 ' 


A>0, 


C>0, 


H2- 


/T 

v « ■ 


A> 0, 


C>o, 


H,^ 


[a 

V 3 ' 


A>0, 


C<0. 



(40) 



Both the solutions with H(t) = Hi = H^ could describe current accerelation era independently on ^, while 
the solution with H{t) = H2 could describe inflationary era independently on A. It has to be stressed that in 
the case ^ > 0, the model effectively describes only one of the two possible phases, because they correspond 
to distinct solutions. 

Now we go back to the general model with arbitrary matter. We first observe that in some cases from the 
generalised Friedmann equation one gets constraints on the possible values of H and p and, as a consequence, 
the corresponding solution is not singular for any finite value of t. In fact, solving (|30p for H^ we obtain 




e>o. 



(41) 



H^ 




4|ei(p + A) 



?<o. 



If <^ > 0, from (|4T|) we get the constraints 



4e 



> p + A> 



h>^H^>0, 



<£.H^ <l, 



C>0. 



(42) 



(43) 



If ^ < equation (|42|) does not give any constraint on H and p. 

Since we would like to describe possible changes of phase during the expansion, for example the exit from 
inflation or the beginning of actual acceleration, we search for the stationary points of the quantity a{t). By 
solving the system of equations (pO| , ([33)1 for H we get 



H' 



(l + w) 



1 



■P 



2(1 - 2^H^) '^ 2(1 - 2ei?2) 

and so a(t) vanishes when H^{t) reaches the (positive) values 



[(1 - 3w)(H'^ + (1 + 3w)H^ - (1 + w)A] 



(44) 



H', 



'fA, 




e>o, 


w = 3, 


1+310 


1 1 /1 1 4eA(l+t«)(l-3tt.) 


w^l/3, 


2{(l-3«;) [ 
l+3w 


1 _L /-I 4|4|A(l+«,)(l-3«.) 


w ^ 1/3. 


[ 2|5|(l-3u.) 



(45) 



This means that when iJ^(to) = Hq > the system goes from a decelerated to an accelerated expansion (or 
viceversa), and this happens only if Ha is in the permitted range of H. If Hq = 0, then there is a singularity 
(a(to) = 00) or a bounce (a(io) = 0, when to < cxi). 



In the case ^ > we must take into account of (P5|) and so, in order to have a possible change of phase, 
the free parameters have to satisfy the fohowing constraints (we consider matter or radiation contributes). 



(A>0,e>0,0,w;=i), 



(A>0,^>0,xi; = 0), 



Hl = lK 



^o = ^(-i + VT+4p;) 



< I CA < 1 , 
1<|?A<2, 

< I CA < 1 , 
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(46) 



l<|eA< 



When A = 0, the only stationary point of d(t) corresponds to Hq = and so the model with (A = 0, ^ > 0) 
can describe only one phase of the expansion. 

In the case ^ < 0, the algebraic equation (|^ does not give any restiction on the values of H^, but 
nevertheless, as we shall see in the following, the parameters ^ and A are not completely arbitrary. Moreover, 
from the third equation in (j45|) we also see that in this case, even for A = there is a non vanishing stationary 
point at i/2 ^ (1 + 3u;)/|^|(l - 3w) with w ^ 1/3. 

For the special case ^ = one can find exact solutions for a{t). In particular, the solutions of ([5n)) - (|5T|) 
with jO > read 



a{t) = ao 



sinh 



1 



3At 



3(1 + ™) 



(47) 



H{t) = coth 



3Ai 



(48) 



p{t) = A sinh 



/3At 



(49) 



ao being an arbitrary constant. As one can see, a{t) vanishes for a finite value of i (i = with the chosen 
initial conditions) and as a consequence both H{t) and p{t) arc divergent at that time. 

As we have seen above, in the absence of matter H{t) — A/3 and a{t) ~ ap exp(-y/A/3t) (de Sitter 
solution) . 

If ^ 7^ we are not able to get explicit solutions for a{t), but nevertheless we can get implicit solutions for 
Hit), which permits to understand the behaviour of the system. 

We start with the two cases (A = 0,^ > 0) and (A = 0,^ < 0). Apart arbitrary integration constants we 
have 



3(1 + w) 


2 


3(1 + w) 



t^l + ^log 
H 2 ^ 



i + V^H 



l-V^H 



t=;^-V^arctan(v^iJ) 



?>0, 



e<o. 



(50) 



(51) 



It has to be noted that in order to have p > 0, in ([SO)) H{t) has to be restricted to the values (,H^ < 1 
in agreement with (|43|) . This means that a{t) does not vanish and the density does not diverge. Equation 
([50)1 effectively corresponds to two different solutions related to the distinct algebraic equations in (HT|) . The 
solution with 1/2 > £^H^ > has an asymptotic behaviour of the kind H{t) ^ 1/t — >■ 0, while the other with 
1/2 < (_H^ < 1 goes as H{t) ^ l/\/C = Const, giving rise to a de Sitter asymptotic behaviour for a{t), that 
is 



a{t) 



floe 



t/Vi 



A = 0. 



C>0, 



2 -^ 



< 1, 



(52) 



As we already said above, the model with such values of the free parameters describes only one phase of the 
expansion. It could be used in order to describe the inflationary era as in Ref. [6[, but it does not provide a 
natural exit from that phase. 



In the second case, equation ([?T|) . H{t) has no restrictions and so there is a singularity for i — >■ (a(i) — >■ 
and p(t) — > oo). This model is not able to describe inflation, because it does not have a rapidly expanding 
phase, but it could describe current acceleration. In fact, putting for simplicity w — in (|45p . we obtain 
^a — Vl^l ^^d using (|44p we see that d{t) is negative or positive according to whether H{t) is smaller or 
greater that Hq- 

Finally we consider the two general cases with (A > 0,i^ > 0) and (A > 0,^ < 0). Apart arbitrary 
integration constants, the implicit solutions for H(t) are given by 



3(1 + U.) 



log 



1 + iJ/^ 



1-iJ/^ 



log 



1 + iJ/^ 



X-HjJW- 



A >0,C> 0, 



where 



a± 



1± 




CA< 



and 



3(1 + w) 



2V^ 



loe 



l + 7?/v/^ 



l-7?/y^ 



/^ 



arctan 



H 



A>0,C<0, 



where 



/3d 




(53) 



(54) 



(55) 



(56) 



Also in this case in order to have p > 0, Hit) in ([5^ has to be restricted to the values a_ < 7?^ < a+ and 
so a(i) does not vanish for finite values of time. The equation describes two distinct solutions, both of them 
giving rise to a de Sitter asymptotic behaviour for a(t) of the kind 



a{i) 
ait) 



aoe- 



floe^ 



/OiJ-t 



lim H{t) = ^a-^ . 
lim H(t) = y/aZ 



(57) 
(58) 



Looking at (|45|) we see that in both the cases w = 1/3 and w = the stationary point is out of the permitted 
range and so also this generalised solution can describe only one phase of the expansion. 

In the last case, equation ([55)1 . p is non negative for H^ > P- and so, in contrast with the analog case 
with A = 0, there is a the solution for which a{t) is always finite for finite values of time. From (|55l) in fact 
we have 



P-<H\t)<(i+. 



^0 < ^ < oo . 



(59) 



13+ <H^{t) < oo. 



oo > t > to . 



(60) 



Both the latter solutions give rise to a de Sitter asymptotic behaviour for a{t) similar to the one in (|57p . but 
with a+ replaced by /3+. 

Depending on the parameters, the previous solutions can describe decelerated and accelerated expansion 
phases. First choosing w = 1/3, from (|^5|) we get 



° 3 



/3_<ij2</3+ if|e|A<| 
H§ > p+ if |e|A > I 



while in the case w = we obtain 



1 



H'o = ^± = 2^ (1 ± v/l-4|e|A 



/3_<i/|</3+ if|^|A<i, 
Hj. > /3+ never . 



(61) 



(62) 
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Using the latter results now we are able to study the changes of phase of the system. For example, if 
w — 1/3 from ^^ it follows that a is negative or positive according to whether H^{t) is higher or lower than 
i^Q = 2A/3. This means that the system passes from a decelerated to an accelerated phase if |^|A < 9/4, 
while it passes from an accelerated to a decelerated phase if |^|A > 9/4. Of course such last situation has to 
be rejected for physical reasons, because both ^ and A have to be small quantiies. 

In the case oi w — 0, the acceleration is negative for i/_ < H{t) < H+ and positive otherwise and this 
means that the model describes three phases, that is 

1) accelerated expansion for /3_ < H^{t) < H^; 

2) decelerated expansion for H^ < H^ (t) < H\ ; 

3) accelerated expansion for H\ < H^{t) < /?+. 

VI. CONCLUDING REMARKS 

In this paper, a toy model of Einstein gravity with a Gauss-Bonnet classically "entropic" term mimicking a 
quantum correction has been investigated. The static black hole solution due to Tomozawa has been recovered 
and generalized with the inclusion of non trivial horizon topology, and its entropy has been evaluated deriving 
the first law from equations of motion. As a result the Bekenstein-Hawking area law has acquired a corrected 
logarithmic area term. A Misner-Sharp expression for the mass of black hole has been found. The same 
model has been used in order to derive Friedmann equations with corrected terms, which reproduce the first 
law with the same formal entropy and energy of the static case, but now related to a general FLRW space- 
time, interpreted as a dynamical cosmological black hole, with a "temperature " associated with Hayward's 
dynamical surface gravity: this result is in agreement with one obtained by the tunneling method applied to 
FLRW space-time in Ref. [33|] . A detailed analysis of all possible cosmological solutions including de Sitter 
and non singular solutions have been provided. These solutions describe several phases, the more interesting 
one is present for ^ < and A > 0. In this case, when w = 0, there is no singularity at i = 0, but an initial 
power-like acceleration, followed by a deceleration phase, and by a final exponential acceleration. 

A generalization of this work may be done by using s-wave approximation and reduction to effective 2nd 
theory as in Ref. 139 1. I t maybe also interesting to generalize the results of this work for F(i?)-gravity (for 
review see Refs. [40,|4l|, where number of BH solutions exist, as in the recent review ^42] ) . 
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